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EMISSION MEASURE DISTRIBUTION IN LOOPS IMPULSIVELY HEATED AT THE FOOT-POINTS 
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ABSTRACT 

This work is prompted by the evidence of sharply peaked emission measure distributions in active stars, and 
by the claims of isothermal loops in solar coronal observations, at variance with the predictions of hydrostatic 
loop models with constant cross-section and uniform heating. We address the problem with loops heated at 
the foot-points. Since steady heating does not allow static loop models solutions, we explore whether pulse- 
heated loops can exist and appear as steady loops, on a time average. We simulate pulse-heated loops, using 
the Palermo-Harvard 1-D hydrodynamic code, for different initial conditions corresponding to typical coronal 
temperatures of stars ranging from intermediate to active {T ^ 3-10 x 10'' K). We find long-lived quasi-steady 
solutions even for heating concentrated at the foot-points over a spatial region of the order of ^ 1/5 of the loop 
half length and broader. These solutions yield an emission measure distribution with a peak at high temperature, 
and the cool side of the peak is as steep as ^ T^, in contrast to the usual ^ T^l^ of hydrostatic models with 
constant cross-section and uniform heating. Such peaks are similar to those found in the emission measure 
distribution of active stars around 10^ K. 

Subject headings: Stars: coronae — emission measure distribution — Sun: coronal loops — X-rays: stars — 
plasmas — Plasma: hydrodynamic modeling 
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1. EVITRODUCTION 

Coronae are important in the study of solar and stellar 
physics for several reasons: they are good tracers of stel- 
lar activity and of dynamo phenomena; also, bright coronae 
identify young stars and stellar formation regions more easily 
than many other characteristics. Nevertheless, fundamental 
aspects of stellar coronae — namely the heating mechanisms 
that sustain the hot, confined plasma — are not well under- 
stood. Although we are not able to observe directly the heat- 
ing processes at work in stellar coronae, the characteristics of 
the observed coronal structures can provide us with an indi- 
rect probe for the properties of the coronal heating. 

The earliest spatially resolved observations of the solar 
corona (e.g. Vaianaetal. 1973) have already showed that 
the hot plasma is highly structured and confined by the 
magnetic field in loop structures, which are considered the 
basic building blocks of the coronae. Considerable ef- 
fort has been devoted in the last three decades to under- 
standing the physics of these structures of confined plasma 
and to develop adequate models that account for the ob- 
served properti es of coronal emission. The fi rst loop mod- 
els (e .g., Rosn eret alJll978i hereafter RTV; IVeseckv et al.l 
ll979trSerio etalJl981[ hereafter S8 1 ) were flux tubes of con- 
stant cross section filled with plasma in hydrostatic equilib- 
rium, and in energy balance under the effects of steady heat- 
ing, heat flux and radiative losses, with the magnetic field 
only confining the plasma. These models demonstrated a 
wide range of validity and satisfactorily reproduced a large 
number of coronal X-ray and EUV observations; both so- 
lar and stellar (e.g. Rosner & Vaiana 1977; Pallavicini et al. 
[l981; Giampapa et al. 1985; Landini et al. 1985; Peres et al. 
[1987; Reale et al. 1988; Reale 2002; Testa et al. 2002). How- 
ever TRACE and SoHO have observed loops apparently in- 
compatible with hydrostatic equilibrium (in terms of their 
spatial distribution of temperature and density) though ap- 
pearing as quasi-static (e.g. Brekke et al. 1997; Warren et al. 
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l2002t iWinebarger et~ani2()0l lGolubll2()0l . Some authors 
claimed that these loops are heat ed non-un iformly and in 
particular at their foot-p oints ( Aschwanden et al...200Q. 20011 
Winebarger et alJl2003l) . On the other hand, if the heating is 
too localized at the foot-points, static loops should be ther- 
mally unstable, as widely discussed sinc e the first wo rks on 
modeli ng of loop structures, e.g. RTV, Antiochos ( 1979), S81, 
iPeres et al. ( 1982) . Recently, detailed analyses of solar obser- 
vations have brought u p again the question of the location of 
the heating release (e.g. lPriest et alJ2000HAschwandeDl2001t 
Reale 2002). 

In order to attempt to explain structures that apparently per- 
sist for time scales longer than the characteristic cooling time 
with foot-point heating models, the models must be dynamic, 
since static solutions with foot-point heating can be unsta- 
ble. Dynamic models of this kind have been recently used 
in several works (e.g. Warren et al. 2002, 2003; Scadaro et alj 
2003; Miiller et al. 2004). For instance. Warren et al. (200j^ 
successfully reproduced several observed characteristics with 
a multi-threaded model, impulsively heating loops to sev- 
eral million degrees and allowing them to cool to TRACE- 
observable temperatures (~ 10* K). 

The study of the stellar coronae allows us to investigate 
the effect of stellar parameters (effective temperature, surface 
gravity, rotation, chemical composition, etc) on the coronal 
models, developed for the Sun. On one hand, when studying 
stellar coronae a first approach is based on the hypothesis that 
the solar corona is an adequate paradigm to interpret the ob- 
servations of stellar coronae. On the other hand, the assump- 
tion of the solar analogy must be validated by comparing the 
characteristics of solar and stellar coronae. 

Recent high quality spectral observations have provided us 
with detailed information on the properties of stellar coronal 
emission. High resolution spectra obtained with, e.g., EUVE, 
in the EUV range, and with the new observatories Chandra 
and XMM-Newton in the X-ray band, allow us to diagnose the 
plasma conditions in a large sample of stellar coronae at dif- 
ferent activity levels. Since stellar coronae cannot be spatially 
resolved by present-day telescopes, we must resort to indirect 
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means in order to compare the properties of coronal structures 
among different stars and with the observed characteristics of 
the solar corona. The emission measure distribution vs. tem- 
perature, EM(T), defined by: 

EM(Ti)= [ nl{T)dV, (1) 
J at; 

where n,, is the electron density, and V is the volume filled 
with plasma at T e AT), contains substantial information 
about the coronal emitting plasma, mainly on its thermal 
structuring, and proves to be a useful diagnostic tool for com- 
paring the coronal emission of different stars. 

The global EM{T) of the X-ray Sun is typically peaked 
at r 2-3 X lO*" K and the ascending cool part is charac- 
terized by a rise as ~ T^^^ i Orlando et alJl2000t iPeres et all 
llpOO), even though its specific properties can vary for differ- 
ent coronal regions and in different phases of activity. Pre- 
vious studies of the solar at mospher e, mostly done in the 
UV band (Ijordan 1980; Brosi us et aljfr9 96: Landi & Landii^ 
led to analogous results. The observed dependence on 
temperature is w ell explained in terms of hydrostatic loops 
iPeres et al .'2001): for a standard RTV uniformly heated loop 
model, the T and n structuring of the plasma along the loop 
yields EM(T) increasing approximately as T^^^. Thus, for 
a corona of optically thin plasma, mostly confined in hydro- 
static loops, this would also yield EM(T) oc T^^^ in the as- 
cending part. The descending part of EM(T) gives us infor- 
mation on the distribution of the hydrostatic loops. Therefore, 
in the stellar case, the analysis of the global EM(T) can pro- 
vide us with information on the structuring of the observed 
corona, and can allow us to test the hypothesis of a corona 
formed by static loops at different maximum temperature. 

Analyses of EM(T) derived from EUV and X-ray spectra 
of several stars have appeared in the recent literature (e.g., 
ISanz-Forcada et ai]l2002l 12003 1 present extensive EM(T) re- 
constructions from EUVE spectra of active stars), and most 
of them show similar features. The EM(T) of active stars 
is typically characterized by an ascending part as r " with 
a > 3/2, and up to ^ 5 in some cases (see e .g. Dupr ee et alJ 
]^93VGriffiths & Jordan'1998"; 'Drake et al."2000; Mewe et alj 
2^01; Sanz-Forcada et al. 2002, 2003; Argiroffi et al. 2003; 
Scelsi et al. 2004). Moreover, the EM{T) of several active 
coronae shows prominent "bum ps", i.e. large amounts of al- 
most isothermal mat erial ("e.^. | Drakelll996|; IPupree 1996, 
I2002h . For example, iDupree etalJ J 199 3') used EUVE data 
for the giant Capella and derived an EM{T) with a well de- 
fined and narrow peak at log(r) '-^ 6.8, and a rise much steeper 
than r^/^; the ana l ysis of X-ray Chand ra spectra of Capella 
dMewe et alJl200UlArgiroffi et alJl2003^ has confirmed these 
findings. Although these results are quite controversial and 
are based on atomic physics parameters still under refine- 
ment, they point to an almost isothermal hot part of the corona 
which, if formed by loops, does not seem to be compatible 
with the predictions of simple static loop models with uni- 
form cross-section. 

Isothermal loops would have fundamental implications for 
the thermodynamics and for the heating of the magnetized 
plasma. In uniformly heated loops, the energy radiated by the 
plasma near to the foot-points (that is at T 10^) is balanced 
by the local heat deposition and by the heat flux from the hot- 
ter plasma further up the loop. However, for foot-point heat- 
ing, the local energy deposition can supply most of the energy 
radiated by the cooler plasma, so that a lower heat flux (shal- 
lower temperature gradient) is needed from the hotter plasma 



and the T profile becomes flatter Since standard hydrostatic 
loop models do not allow stable solutions for heating deposi- 
tion over a spatial region smaller than about L/3 (e.g. S81), 
we decided to explore the possibility of obtaining stable loops 
with a non-constant and thus dynamic heating. 

The possibility of a corona composed of loops heated 
episodically, close to their foot-points, poses two basic ques- 
tions: 

• Under what conditions can such loops be stable for time 
scales much longer than the characteristic cooling time? 

• If long-lived, foot-point heated loops exist, is their 
EM(T) different from that of uniformly heated loops? 

In this work we model the properties of foot-point heated 
coronal loops, using a 1-D hydrodynamic model. We 
run simulations with pulse-heating concentrated at the foot- 
points, and investigate the existence of long-lived solu- 
tions; we study the EM{T) associated with these solu- 
tions. In particular, we investigate the existence of sta- 
ble loops characterized by an EM(T) with a narrow peak 
and a slope steeper than T^^^, i.e. characteristic of those 
derived from observatio ns of active stars . We note that, 
for the solar case, CargiU & KlimchukI (12004, see also 
Cargill & KHmchuk 1997; KHmchuk & Cargifl 2001) demon- 
strated that the EM(T) provides a useful diagnostic for peri- 
odically heated plasma, by using a zero-dimensional model 
(i.e., each loop is characterized by a single, averaged, value 
of density and temperature) that considers nanoflare heating 
(Cargill 1994). 

In ^we describe the loop modeling and the properties of 
the simulations. In ^we analyze the results of the simula- 
tions. In l^we discuss our results and their implications, and 
then we draw conclusions from our work. 

2. LOOP MODEL AND SET OF SIMULATIONS 

We model the plasma confined in a single magnetic flux 
tube. We will discuss later how this single loop model can be 
useful for interpreting the overall emission from a corona. 

In our model, the magnetic field has only the role of con- 
fining the plasma; thermal conduction is effective only along 
the magnetic field lines. The model assumes constant cross- 
section along the loop. In our study, the loop half length is 
assumed to be L = 10"^ cm, as observed for relatively long 
loops on the Sun. We discuss below the implications of our 
assumptions. 

We simulate coronal loops heated close to their foot-points 
by episodic pulses. In order to investigate the effect of some 
basic parameters on the solutions for the loop plasma, we have 
reduced the number of free parameters as much as possible. 
For instance, we assume identical conditions in both loop legs 
and, therefore the loop model is symmetric about the loop 
apex. Because of this symmetry, the equations for the loop 
plasma are solved for half of the loop only. 

The simulations are long lasting to search for long-term sta- 
bility, i.e. for loops settling in a state steady on average, if such 
solutions exist. 

We use the Palermo-Harvard code ( iPeres et al] Il982t 
Betta et al. 1997); a 1-D hydrodynamic code that consistently 
solves the time-dependent density, momentum and energy 
equations for the plasma confined by the magnetic field: 
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TABLE 1 
Parameters of the simulations. 



model " 






Initial conditions 




Heat pulses 








7^ d 

-* max 


Phase " Eq ' 




{Eh)' 


'run 


C 


1 


3 


1 0.45 


--2200 


L/3, L/5, L/10 £o, 4£o Tcooi/4, Tcooi/2 


10000 


H 


1 


10 


36 30 


~ 1200 


L/3, L/5, L/10 £o, 4£o r^ooll^, r^ooi/^ 


5 000 



and H indicate the models of cooler (T = 3 X 10* K) and hotter loops (T = lO' K) respectively. 
''The set of simulations considers all possible combinations of the values of the parameters (a, {Eh), t). 
'' Loop semilength in units of lO'" cm. 

''Maximum temperature (= apex temperature) in units of 10* K. 
"Base pressure in units of dyn/cm^. 

'Heating per unit time and per unit volume (erg cm"^^ s"') from RTV loop scaling laws. 
^ Loop cooling time (see Eq.|6| in seconds. 

''Gaussian parameter of the spatial extent of the heating (see Eg. 1101 . 

'Heating intensity averaged in time and along the loop, expressed in terms of the static heating, Eq. 
'Period of the heat pulses. 
■"Total simulation time (s). 



nniii 



dv dp 



dt 



(3) 



dt OS as OS OS 

with p and e defined by: 



p = {\ + (3)nK^T 



e=-p + nPx, 



(5) 



where n is the hydrogen number density; s is the spatial co- 
ordinate along the loop; v is the plasma velocity; niu is the 
mass of hydrogen atom; /i is the effective plasma viscosity; 
PiT) represents the radiative losses function per unit emis- 
sion measure; f3 is the fractio nal ioniz ation, i.e. rif^/riu', k is 
the Spitzer conductivity (Spit zeill962h : is the Boltzmann 
constant; and x is the hydrogen ionization potential. Eu is 
an ad hoc heating function of both space and time; this is the 
main parameter we vary to study the characteristics of the so- 
lutions, and it will be described in detail in ^2.21 The numeri- 
cal code uses an adaptive spatial grid to follow adequately the 
evolving profiles of the physical quantities, which can vary 
dramatically in the transition region. 

2.1. Initial Conditions 

As initial conditions, we consider hydrostatic loop solutions 
of the S81 model with uniform heating. In particular we have 
selected two different solutions with maximum temperatures: 
Tinax = 3 10'' K and Tmax = 10^ K. Th e initial model atmo- 
sphere uses the IVernazza et alJ Jl976h model to extend the 
S81 static model to chromospheric temperatures (the mini- 
mum temperature is Tmin = 4.4 • lO"* K). Table ^ summarizes 
the characteristics of the initial loop conditions: maximum 
temperature, Tmax; base pressure, phase', heating per unit time 
and per unit volume, Eq; and the characteristic cooling time, 
Tcooi, according to ^rio et al. ( 1991): 



Tcool ~ 120' 



(6) 



This parameter is important in many respects, e.g. for com- 
parison with the time interval between heat pulses, and other 
parameters of the simulations as discussed below. 



Note that the parameters of the loop in Table [J satisfy the 
scaling laws derived from the static RTV loop model: 



Pbase 



1.4-103 



1 
L 



Eo ^ W 



„I.I7i--().83 
■ Fbase^ 



(7) 



(8) 



Our aim is to find solutions corresponding to loops in steady 
state conditions over long time scales, as observed. We find 
that, in general, the initial conditions have little influence on 
the loop evolution driven by the repeated impulsive heating, 
and that they can be important only in a region of the parame- 
ter space at the boundary between stability and instability. In 
particular, the results can change if the loop is initially already 
hot and dense. 

2.2. The Heating Function 

The heating function, Eii{s,t), is assumed to be a separate 
function of space and time: 

En(s,t) = Ho-g(s)-f(t) (9) 

Amount of Energy Release One of our main goals is 

to model stellar coronae of activity levels from intermediate 
to high, showing features not reproduced by standard mod- 
els. The characteristic plasma temperatures for such coronae 
range from a few million degrees up to > 10^ K. 

The reference value for the time-averaged intensity of the 
heating is the energy required to heat the initial hydrostatic 
model (£■() of Table^. We run two sets of simulations choos- 
ing Ho such that the heating, averaged in time and along 
the loop, (Eli) (= ^7 Jy JfEu(s,t)dVdt), is equal to £0 corre- 
sponding to Trtv = 3 X 10^ K and Trtv = 10^ K respectively^. 
We also consider higher values of the heating i{En) = 4£'o) 
because: 

1 . we want to check whether unstable loops become stable 
with stronger heating; 

2. foot-point heated loops are cooler than uniformly 
heated loops with the same total energy input; we in- 
clude simulations with stronger heating to compensate 
for this effect. 

^ Trtv corresponds to Tmm of E^^ 0. Trtv and Eo are linked through the 
scahng laws of Equations Q and js). 
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Spatial Distribution of Heating The spatial distribution of 

the heating is described by a Gaussian function: 

g(i) = e-(--^o)V(2<T=)^ (10) 

centered at the loop foot-points (i.e. sq = 0) for all the simula- 
tions. 

We assume that the heat deposition is exactly the same in 
both legs of the loop, thus, as mentioned above, we have sym- 
metry about the loop apex. 

The spatial extent of the heating is a fundamental param- 
eter of our study, given our interest in the effect of the con- 
centration of the heating on the stability and on other charac- 
teristics of the solutions. A preliminary set of runs has indi- 
cated that the boundary between stable and unstable solutions 
is cr ~ L/5. Therefore, we will discuss in detail the results for 
a = L/10,L/5,L/3, i.e. smaller than, equal to and larger than 
the critical value. For a > L/3 the solutions do not depart 
significantly from the static solutions. 

Temporal Distribution of Heating Pulses Given the tran- 
sient nature of the simulated heating, a natural reference point 
is represented by the loop models with heating that produces 
flare events. In these models the heating typically consists of 
two terms: a constant and uniform heating keeping the steady 
conditions of the initial loop, i.e. the steady coronal heating; 
and a transient much larger heating which triggers the flare. 
In our analysis we neglect the constant term and assume no 
steady coronal heating. 

The heat pulses are periodic with period r, and duty cycle 
10%, i.e. the heat is active only for a pulse lasting 1/10 of 
the period. Each heat pulse is a step function, constant when 
active and zero otherwise. We consider values of r smaller 
than the cooling time of the loop, TcooI, in order to prevent the 
loop from collapsing, but long enough to induce significant 
changes in the loop. 

Our choice of periodic pulses, instead of a random distri- 
bution in time, allows us to limit the free parameters and to 
focus on the most critical ones. We have also run test cases 
with rand om heat pulses in o rder to check for possible dif- 
ferences ( Testa & Peresl2003h : solutions do not differ signifi- 
cantly for similar pulse parameter values (i.e. ave rage energy 
release and average interval between pulses). Also lPeres et al.l 
{1993) have shown that the results do not differ substantially 
from those obtained with random pulses with equal average 
repetition time and duty cycle. 

The simulations are run for at least ~ 10 periods and a total 
time, fiun, much longer than the cooling time, in order to be 
able to distinguish between stable and unstable solutions. The 
total times are of the order of several hours, comparable with 
the time coverage of actual coronal observations. 

We summarize the parameters of the simulations in Tabled 

3. RESULTS 
3.1. Evolution of Solutions 

We know from previous works that loops continuously 
heated at the foot-points are not stable, for a < (Tcriticai- For 
comparison with pulse-heated loops, we first show the results 
for a loop with continuous heating. 

Figure [3 shows the evolution of temperature and density 
for a loop evolving from initial hydrostatic conditions of 
Tmax ^ 10^ K, with a foot-point heating {a = L/3) constant 
in time. The evolution is shown as a 3D plot, where the tem- 
perature, T, and the density, n, are plotted vs. time and vs. the 
coordinate, s, along the loop, s = 0, corresponds to the loop 



base, s = 10'" cm, to the apex; the 3D plot is oriented with 
the initial profiles toward the observer. The figure shows that 
the plasma progressively cools down and condensates at the 
loop apex. As discussed by e.g. RTV and S81, a configuration 
with temperature and density inversion {s{Tjnax), '?(nmin) < L) 
is thermally unstable. Despite the continuous energy release, 
the loop collapses on time scales only slightly longer than the 
characteristic cooling time. 

Figure 121 and |3l show results for pulsed heating, and in par- 
ticular for an unstable and a stable loop respectively, in the 
same format of Figure^ 

The solution shown in Figure |2l is for a loop with semi- 
length L = 10'*^ cm and heating function parameters: r = 
^0001/2; cr = L/IO; and (Eu) = 4£'o (see ^or Table[l]for their 
definition). Initially, the solution does not depart markedly 
from average conditions, but after t ^ 2T(;oo1 the instability 
develops at the loop apex, i.e. the least heated place in the 
loop. As the apex plasma cools slightly, the radiative losses 
increase. Since the apex plasma is not sustained by significant 
local heat deposition, it cools down even more and condenses; 
the increase in density and the decrease in temperature fur- 
ther enhances the radiative losses. Once the instability is trig- 
gered, the loop quickly collapses because the radiative losses 
increase for decreasing temperature and increasing density. 

All the unstable cases show similar temperature and density 
evolution; the more concentrated the heating, the faster the 
instabilityoccurs. The instability invariably starts at the apex. 

Figure |3] shows an example of a long-lived solution; for a 
more direct comparison we show the loop with the same char- 
acteristics as that shown in Figure|2l except that the spatial de- 
position of the heating is cr = L/5 instead of L/10 (Figure|2jl. 
In spite of the fluctuations due to the heat pulses, the loop 
is stable on long time scales and settles to a state of higher 
density because of the significant chromospheric evaporation 
driven by the heat pulses. 

Figure |3 shows the evolution of the plasma temperature 
and density at the loop apex for the loop models with pulsed 
heating as compared with the solution for constant heating. 
The evolution of the plasma properties at the loop apex are 
presented for the loop constantly heated at the foot-point 
(cr = L/3; dashed lines), and for the loops subject to heat 
pulses with different spatial distribution: a = L/3 (thick light 
grey); a = L/5 (solid black); and a = L/10 (dark grey). The 
figure shows the instability of the constantly heated loop and 
of the solution with heat too concentrated at the foot-points 
(cr = L/10). The two solutions with a = L/5 and a = L/3 ap- 
pear more stable. 

Effect of Heating Parameters The critical parameter for 

loop stability is the spatial width of the heating function: we 
have stable solutions for cr > L/3 and unstable solutions for 
cr < L/10; for a ^ L/5 the loops are on the edge of stability, 
and other characteristics of the heating (intensity, interval be- 
tween pulses) become important. We find consistent results 
for both the cooler and the hotter solutions. 

The interval between pulses, t, does not appear to be a 
critical parameter for stability as long as it is not too small 
(t <C Tcooi), i.e. too close to the constant heating case that is 
unstable for cr = L/3 (see Fig.Q, nor too long with respect to 
the cooling time, since the loop would then catastrophically 
cool. 

3.2. Emission Measure Distribution 
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Fig. 1. — Evolution of the temperature (left) and the density (right) distributions along (half of) a loop with semilength L = 10 cm, initial maximum 
temperature Tmax = lO' K, and constantly heated at its foot-points with a = L/3. T (in K) and n (in lO' cm"^) are plotted vs. time (in seconds) and vs. the 
coordinate along the loop (in units of lO' cm). The loop top is on the left, and the 3D plots are oriented with the initial profiles toward the observer. 




Fig. 2. — Evolution of the temperature (left) and the density (right) distributions along (half of) an unstable pulse-heated loop. The heating function parameters 
are: r = Tcooi/2, u = L/ 10, (£h) = 4£o (see TablefTlfor their definition). The 3D plots are in the same format as in FigurefTI 



Besides investigating the stability of the solutions, our mod- 
eling effort aims at studying the effect of shrinking the heating 
region at the loop foot-points on the EM{T), which is one of 
the main derived quantities that allow us to study and to com- 
pare the coronal properties of different stars. 

As described in section |2j we analyze two different sets of 
models with maximum temperatures of about ^ 3 x 10^ K, 
and ^ 10^ K, respectively. Since we find results that are qual- 
itatively similar in both cases, we will discuss in detail the 
model of hotter loop, because one of the questions we want 
to address is the possibility of reproducing the EM(T) of ac- 
tive stars (generally characterized by peak temperature of the 



order of 10^ K) with the proposed model. 

From our simulations we can model the global emission of 
a stellar corona composed of many loops impulsively heated 
by microflares occurring close to their foot-points. Under the 
assumption that all of the loops have a statistically analogous 
evolution, taking a sample of profiles at different times in the 
evolution of a single loop, as shown for example in Fig.|3] is 
equivalent to observing all of the loops simultaneously, with 
each loop at a different stage in its evolution. We take 200 out- 
puts from each simulation, uniformly sampled in time, and 
consider each output as a snapshot of an independent loop. 
Thus, the total emission from the corona, composed of all 
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Fig. 3. — Evolution of the temperature (left) and the density (right) distributions for a stable pulse-heated loop. The heating function parameters are: t = T;.ooi/2, 
a = L/5, {Eii) = 4Eq, i.e. the parameters are the same as those of the loop of Figurel^ except for cr. The 3D plots are in the same format as in FigurelTI 
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Fig. 4. — Evolution of the temperature (Teft)^ana the density (right) at the loop top, for different spatial and temporal characteristics of heating: continuous 
heating concentrated at the foot-point with spatial scale cr = L/3 (dashed lines), periodic heat pulses with cr = L/3 (thick light grey), cr = L/5 (solid black), 
cr = L/10 (dark grey). 



these loops, is the sum of the emission of all the snapshots. 

From the temperature and density profiles of each out- 
put we derive EM{T) and then we sum them to obtain the 
global EM{T). We derive EM(T) in each temperature bin. 
At;, as EM(Ti) = J2j'''eU) dsj where j spans over the spatial 
bins whose temperature falls within AT}, and the temperature 
bins are constant in logT such that AlogT = 0.1. Therefore 
EM(T) obtained with this procedure differs from the defini- 
tion of Eq.([0 by a normalization factor that corresponds to 
the cross-sectional area of the loop. 

Figure |5] shows EM(T) for the loop models with average 
coronal temperatures ~ 10^ K, and heating functions with 



different spatial widths. The solutions with cr = L/3 (top), 
cr = L/5 (middle), and ct = L/10 (bottom), show the effect on 
EM(T) of varying the heating concentration. The lines in the 
lower part of the plots of Figure |5] correspond to EM(T) for 
each output represented with lines darkening from t = (light 
gray) to f = fmn (dark gray). The upper thick black line repre- 
sents the total EM(T), which is the sum of all the individual 
distributions. Therefore the total EM(T) is the one that would 
characterize a corona composed entirely of such loops. In or- 
der to allow for an easier comparison of the static model with 
the dynamic simulations, EM(T) of the initial static solution 
(dashed-dotted line) is also shown arbitrarily shifted next to 
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the total EM(T). In Figure |5]the power laws corresponding 

c>c are also plotted as useful 
points of reference, as discussed in ^ 
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5. — Emission measure distribution, EM(T), for tlie solutions with 
Tmax ~ WK, {Eh) = 4£o, and cr = L/3 (top), a = L/5 (middle), a = L/10 
(bottom). The lines in the lower part of the plots correspond to the EM(T) of 
the individual outputs, regularly sampling the whole loop evolution; the lines 
get darker from ? = (lighter gray) to ? = ;,u„ (darker gray Hnes). The dashed- 
dotted lines correspond to the EM(T) of the initial static condition, while the 
upper thick black line is the total EM(T). For comparison the dotted lines 
indicate the power laws coiresponding to EM(T) oc T'^l'^ and EM(T) oc . 



The EMiT) of the initial hydrostatic solution is quite flat, 
close to the T'^l^ power-law. Figure |5] shows that EM(T) can 
clearly be used to characterize the dynamic simulations since 
we see that the total EMiT) changes consistently as the heat- 
ing becomes more concentrated in the foot-point region. The 
main difference of the total EMiT) (solid lines) from the stan- 
dard hydrostatic distribution is the presence of a well-defined 
peak at high temperatures i^ \if K). This is also the temper- 
ature range corresponding to the bulk of the emission from 
active stellar coronae. 

For (T = L/3 the EMiT) is close to the initial curve. Narrow- 
ing the region of energy release to ct = L/5 causes the EMiT) 
to steepen and approach the scaling derived from stellar ob- 
servations. The more concentrated the heating, the wider the 
peak. 

The EMiT) does not depend on the parameter r as it does 
on the concentration of the heating. 

In Figure|6lwe compare the EMiT) derived from our model, 
for a = L/5, to an EMiT) derived by Scelsi et al. (2004) from 
XMM-Newton spectral observations of the active star 3 1 Com. 
The EMiT) derived from a RTV hydrostatic model is also 
shown. The EMiT) derived from our model, with impulsive 
foot-point heating, is qualitatively similar to the EMiT) of 
31 Com. 

4. DISCUSSION AND CONCLUSIONS 

Coronal loops are known to be unstable if heated con- 
tinuously in a sufficiently small region at their foot-points. 
Also, the EMiT) of stable, uniformly-heated loops of con- 
stant cross-section is known to obey a power law with index 
~ 3/2. In the present work, we have investigated the existence 
of long-lived coronal loops heated by sequences of pulses lo- 
cated at their foot-points and studied the changes induced in 
the EMiT) by this form of heating. In particular, we exam- 
ined relatively long (2 x 10'" cm) and hot (w 3 x 10^ K and 
10^ K) loops appropriate to investigating both the solar envi- 
ronment, which we can spatially resolve and compare directly 
with models, and the coronae of active stars, where evidence 
of different EMiT) has been found. 

We find long-lived dynamically stable solutions in cases 
with heating concentrated close to the foot-points, with val- 
ues of a down to ^ L/5. Therefore we find stable solutions 
for much more concentrated heating compared with the solu- 
tions of the S81 model; the latter, using a heating spatial dis- 
tribution of the form Hoe'"!"^, does not yield stable solutions 
for ^ L/3. The value of ct ^ L/5 represents the border 
zone between rapidly unstable and long-lived loops for both 
sets of simulations with different maximum temperatures (see 
Tab. [0, provided that the heating is strong enough to sustain 
the loop. We note, however, that close to this boundary some 
unstable cases appear quasi-steady when observed over time 
scales that are long compared to the cooling time. We find that 
the time scales of the loops' evolution do not critically depend 
on the intensity of the heating, or on the period of heat pulses, 
with the caveats discussed in the previous section. 

Our work has therefore proven that long-lived footpoint- 
heated loops can exist, provided that the heating is intermit- 
tent with the appropriate period - a fraction of the loop cool- 
ing time - which should be long enough to depart from con- 
tinuous heating, but shorter than the loop plasma cooling time. 
The intermittent heating allows the plasma at the loop apex to 
drain downward to the chromospheric region (as confirmed 
by the downward velocity of the plasma close to the apex 
in between the pulses), and prevents the accumulation of the 
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EM{T) QC (dotted line) and EM{T) oc T^/^ (dash-dotted line) are shown as well. EM{T) from models are arbitrarily scaled. 



plasma at the loop top, and thus the thermal instability. Our 
results are in agreement with those of other works modeling 
footpoint-heated loops (e.g., M iiller et al..,2004) . 

In the present work we also showed that the EM(T) of foot- 
point heated loops are significantly different to those of con- 
ventional hydrostatic loops. The EM{T) of foot-point heated 
loops shows a well-defined peak, which becomes wider and 
wider as the heating becomes more concentrated towards the 
foot-points. This holds true for both sets of simulations at 
different temperatures analyzed here (see Tab. The low 
temperature side of the peaks has a steeper slope than the 
EM(T) cx r^/^ of the static case. The slope found for the 
stable cases with smaller a approaches the power law of 
EM{T) oc T^. We note, in passing, that loops with coro- 
nal cross-section larger than in their chromospheres may give 
similar results (e.g., Schriiver et al. 1989; Ciaravella et al. 
[T996: Sim & Jordan 20031) . We plan to investigate this pos- 
sibility in the future. 

The EM(T) of the impulsively heated loops modeled here 
show a qualitative agreement with those derived from sev- 
eral X-ray and EUV spectroscopic observations of active 
stars. Recent analyses outline a scenario in which hotter {T ~ 
10^ K) plasma is characterized by physical conditions that 
are fundamentally different from the cooler {T <3 x 10^ K) 
plasma, which probably belongs to a different class of struc- 



tures. In particular, there is increasing observational evidence 
of hot plasma with density two orders of magnitude higher 
than the cooler plas ma. Th ese density and filling factor re- 
sults (e.g.. Testa et al,||2n03) support a scenario in which, for 
increasing activity levels, significant flaring activity may be 
present, yielding the hotter plasma structures. In such a sce- 
nario, it is likely that hot loops are sustained by impulsive 
energy release. Therefore, the theoretical model discussed 
here is based on assumptions that are consistent with the ob- 
servational evidence, while reproducing the steep and peaked 
EM{T) widely found for active stars. 

There is room for several improvements to be made to our 
model: we plan to model loops with a larger cross-section 
in the corona; and we will consider the larger set of models 
for unstable, foot-point heated loops, since the superposition 
of several loops of this kind may help to explain unresolved, 
monolithic loop structures composed of many strands, or even 
whole stellar coronae. 
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